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MULTI LAND: a neutral landscape generator 
designed for theoretical studies 

Lionel Roques 


1 Presentation of the software 


The main goal of Multiland is to generate neutral landscapes made of several 
types of regions, with an exact control of the proportions occupied by each type 
of region (see Gaucherel et ah, 2014, for a review of current landscape modelling 
approaches). An important feature of the software is that it allows a control of 
the landscape fragmentation. It is intended to theoretical studies on the effect 
of landscape structure in applied sciences. The user-friendly interface contains 6 
entries (Fig. [I]). 



Figure 1: User interface of the Multiland software. 
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Inputs 


Landscape parameters: 


• grid size: an integer k > 2. The landscape is described by a lattice of size 
k x k; 


• proportions: a vector V = (vi,...,%) of length M > 2. The length M 
of the vector gives the number of different regions in the landscape. The 
proportion of the lattice occupied by the region associated with a value p G 

[1 , M] is equal to p p = v p / ( I whenever k 2 p p is an integer. Otherwise, 


0 = 1 / 


the number of cells n p corresponding to region p < M is given by the floor of 

M—l 

k 2 n p : [ k 2 p p \. The number of cells in the last region is given by k 2 — n p- 

p =i 


• fragmentation index: a real in [0,1]. The value 0 corresponds to the most 
aggregated landscapes; the value 1 corresponds to the most fragmented 
landscapes. 


Computation parameters: 

• max. comput. time: number of seconds allowed for the generation of 1 
landscape. Default value is 5s. Computation time depends on the grid 
size; on the value of the fragmentation index (longer times are necessary for 
values close to 0 or 1); and on the tolerance (see below); 

• tolerance: the computation is stopped when the maximum computation 
time is reached or when the algorithm has generated a landscape with the 
desired fragmentation index ± the tolerance. 


Landscape generation: 

• Number of landscapes: once the button “Generate landscape(s)” is pressed, 
the indicated number of landscapes is generated. The matrixes describing 
the landscapes are save in a file “U.txt”. The corresponding fragmentation 
indexes are save in a file “fr.txt”. 
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Figure 2: The 4-neighbourhood system: an element £ of S and its four neighbours 
i = 1...4. 

2 Theoretical aspects 

The software Multiland is based on a stochastic landscape model with M types 
of regions. This model should: 

- allow an exact control of the proportion attributed to each region; 

- lead to more or less fragmented landscapes, with a control of the degree of 
fragmentation. 


State space. Let S be the lattice made of the of the cells i G {1,..., k}, 

i G (1,..., k}, equipped with a neighboring system VL. We consider the random 
field oo defined on the lattice S. This random field assigns each site £ = (i, j) 
(( 6 E) a value p in {1,2,..., M}, while the number n p = card{ca^ = p, £ G 5}, 
takes a pre-fixed value. Under these assumptions let U be the corresponding state 
space. 

Neighborhood system. We adopt a 4-neighborhood system as shown in Fig. [2} 
In order to take account of the periodicity of the environment, the domain is 
considered wrapped on a torus. 

Fragmentation index. For each oo in Q, we set 

= \ = 
r^es, 

the number of pairs of neighbors £, rj G S such that oo% = oo v (!{•} is the indicator 
function). The statistic of the model, s(oo), is directly linked to the habitat 
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fragmentation: a landscape pattern is all the more aggregated as s((n) is high, 
and all the more fragmented as s(cu) is small. This statistic will be used below in 
the construction of the landscape model. 

In order to obtain a fragmentation index which belongs to [0,1] (0: most 
aggregated landscapes, 1: most fragmented landscapes), we proceed as follows. 
First, we choose q in {1,2,..., M } such that n q > n p for all p G {1,2,..., n). 
The region associated with the value q corresponds to the “background” of our 
landscape. Then, for each other region z associated with a value in / = {1,..., q— 
1, q + 1,... n}, we define the fragmentation level 


s z (u}) 

B(n z y 


( 2 . 1 ) 


where s z (u ) 


l = U 'n 

»?>£eE, ri£V^ 


z}, corresponds to the number of pairs of 


similar neighbors in the region z and B{n z ) = 2n z — ccil(2 Ak /nT) corresponds to 
the maximum possible number of similar neighbors in a region made of n z cells. 
Thus, f z € [0,1]. Finally, we set 


z£l 


( 2 . 2 ) 


and we say that J r (ui) is the fragmentation index of the landscape u. 


Gibbs measure. The landscape model is based on the Gibbs measure P, defined 
over 11 by 

P(X = u>) = ±e~ u M, (2.3) 

with the partition function Z = e ~ (/ G) anc [ ^j ie energy function U : 

VL -> M, 


U(u ) = ~Ps(u), 

with (3 G M; thus, the probability density becomes 


(2.4) 


P(X =u) = (2.5) 

Zj 

The proposed landscape model favours the pairs of neighbours with the same 
value, while /3 increases. Hence, f3 directly controls the topology of the landscape 
patterns. 


Generation of the landscapes. Samples from the distribution (2.5) are ob¬ 
tained by using a Metropolis-Hastings algorithm. More precisely, starting from 
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a random initial landscape, a Metropolis-Hastings algorithm is run until (1) the 
maximum allowed computation time is reached or (2) the distance between the 
fragmentation index of the current state and the desired fragmentation index is 
smaller than the tolerance (see Section 1). The parameter f3 is adjusted during 
the algorithm in order to increase or decrease the fragmentation index, depending 
on the desired fragmentation index. 


3 References and contact 


A simpler version of the algorithm was used in Roques and Stoica (2007) to gen¬ 


erate binary landscapes. This new version has been developed in the framework 
of the PEERLESS ANR project “Predictive Ecological Engineering for Land¬ 
scape Ecosystem Services and Sustainability”. For more information, contact the 
author at lionel.roques@avignon.inra.fr. 


4 Downloading the software 

The Matlab® source code is available at http://multiland.biosp.org. 
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